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Abstract 

We study rooted spiral trees in 2, 3 and 4 dimensions on a hyper cubic lattice 
using exact enumeration and Monte-Carlo techniques. On the square lattice, we also 
obtain exact lower bound of 1.93565 on the growth constant A. Series expansions 
give 9 = -1.3667±0.0010 and v = 0.6574±0.0010 on a square lattice. With Monte- 
Carlo simulations we get the estimates as 9 = — 1.364±0.010, and v = 0.656±0.010. 
These results are numerical evidence against earlier proposed dimensional reduction 
by four in this problem. In dimensions higher than two, the spiral constraint can 
be implemented in two ways. In either case, our series expansion results do not 
support the proposed dimensional reduction. 

Keywords: Spiral trees, Exact enumeration, Dimensional reduction 

Spiral structures are very common in nature. Some examples of the beautiful spiral 
structures in galaxies, shoot arrangement in plants, polymers with spiral structure etc 
may be found in the book by Hargittai [1]. In statistical mechanics, lattice models of 
spiral self avoiding walks have been studied and can be solved exactly in two dimension 
[2, 3], though no solution is known for the self avoiding walks without the spiral constraint. 
A model of spiral trees and animals was proposed by Li and Zhou [4], which based on 
numerical studies was suggested to be in a new universality class. This problem was 
further studied by Bose et. al [5]. Based on the numerical evidence, and guided by the 
fact that magnetic field acting perpendicular to the motion of a charged particle produces 
spiralling motion and reduction by two in effective dimensionality, they conjectured that 
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spiral tree problem would show a dimensional reduction by four. They conjectured the 
exponents of the spiral tree problem follow the following relations 



= (d-A)u pl for d = 2 
9 = (d-4)v ± for d > 2 



(1) 
(2) 



where 6 is the entropic exponent and v v i and v± are the exponents related to the radius 
of gyration in the plane in which the tree has a rotational constraint and perpendicular 
to that plane respectively. 

Since then this problem has not been studied further. Dimensional reduction is an 
intriguing possibility. The lattice tree model without spiral constraint is known to show 
a dimensional reduction by two [6]. The directed version, show a dimensional reduction 
by one. For both models, the tree and animals are believed to lie in the same universality 
class. In this paper, we revisit the problem and obtain a significantly longer series for 
rooted spiral trees. Specifically in two dimensions we have added twelve terms to the 
earlier series of 25 terms. In three and four dimensions, we generated a seventeen and a 
thirteen term series respectively. The earlier known series in three and four dimensions 
had thirteen and nine terms. In the process, we also correct some mistakes in the earlier 
reported series. We also perform Monte-Carlo (MC) simulations using the improved 
incomplete enumeration algorithm [7] and generate spiral trees up to sizes of 1000 in two 
dimensions. Our analysis of exact series and MC samples do not support the conjectured 
dimensional reduction by four in this problem. 

A lattice tree is a cluster of connected sites which contains no loops. Spiral trees are 
a subclass of lattice trees. In a tree every cluster site is attached to the origin through a 
unique path. In a spiral tree, this path has a specific rotational constraint. 

We define rooted spiral tree as a acyclic connected subgraph of a lattice such that the 
projection of the path joining any site of the tree to the root on x — y plane contains 
no left turn (Fig. 1). We will measure the size of a spiral tree by the number of sites 
present in the tree. These are called spiral site trees. The number of possibilities of spiral 
bond trees are more than that for spiral site trees but both are thought to lie in the same 
universality class. For example, the site marked as X in Fig.l is not allowed in spiral site 
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tree as it would introduce a loop. But it can be present in a spiral bond tree. 

Let the total number of distinct rooted spiral trees be A n . This is expected to have a 
asymptotic behaviour of the form 

A n ~ C\ n n- e (3) 

where C is a constant, 9 is a critical exponent and is expected to depend only on the 
dimension of the lattice and A is known as the lattice dependent growth constant. The 
existence of growth constant A for unrooted lattice trees and animals has been proved 
rigorously using concatenation and super multiplicity arguments [8]. Also a rigorous 
lower bound for 9 for unrooted lattice trees and animals has been proved [9] using pattern 
theorem. Specifically, it is 9 > (d — l)/d, for any dimension d > 2. Eq. 3 is expected to 
hold for most cluster enumeration problems on regular lattices, though other asymptotic 
forms are also possible. For example, for spiral self-avoiding walks on square and triangular 
lattices, A n tends to a stretched exponential in n in the asymptotic limit [2, 3]. Though 
we do not prove existence of A and 9 for spiral trees in this paper, we derive a lower bound 
for A. Also, since spiral trees are a subset of lattice trees, X sp irai < A a «, where X a u is the 
growth constant for all trees, \ a ii ~ 3.795 on a square lattice [10]. In two dimensions, 
we have derived the generating function for enumeration of a subset of all possible spiral 
trees. The value of growth constant for this subset is 1.93565. This gives a lower bound 
on the growth constant \ sp irai of the spiral trees on a square lattice. This bound will be 
derived in Section 1.1. 

For the conventional lattice animals, one can prove 9 > through concatenation 
and super-multiplicity arguments [8, 9]. Concatenation does not work for spiral trees. 
Interestingly, our numerical studies give a negative value of 9 in two and three dimensions. 

The spiral trees are anisotropic. We measure the average extent of a n-site spiral tree 
in x — y plane and perpendicular to the x — y plane through the moment of inertia tensors, 
I p i t n and 7_L,n respectively. In the asymptotic limit, they are expected to vary as 

I pl , n ~ An 2 ^ +1 (4) 

and 
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I±, n ~ An 2 ^ +1 (5) 

where u pi and define the length scale of the spiral tree in planar and perpendicular 
direction respectively. 

1 Two dimensional lattice spiral trees 

Some pictures of randomly generated large spiral trees are shown in Fig. 2 (details later). 
One notes very long one dimensional structures with infrequent turns. Hence, simple 
counting of structures of kind shown in Fig. 3 should give a good estimate of the growth 
constant A. The generating function of trees of this type is easy to determine. If Ai(x) is 
the generating function, we get 

which gives Ai(x) = 1+ ^ l J^]_ x 3 ■ The number of trees of this type grows as Ai n , with A x = 
1.754878. It is straightforward to include more complicated branches in such counting to 
get a better lower bound. This we proceed to do below. 

1.1 Lower Bound on Growth Constant on Square lattice 

Consider a subset of all the spiral trees on a square lattice rooted at the origin, which 
lie strictly in the first quadrant x > 0, y > 0; starting at the origin, and not touching 
y = and y — 1 except at points (0,0) and (0, 1) respectively. If Q(x) is the generating 
function for spiral trees in a quadrant and if g 4 n is the coefficient of x n in the expansion 
of {[Q{x)} 4 )/x 3 , then 

A n > g 4 ,« (7) 

where A n is the n th term of A(x), the generating function of all spiral trees on the square 
lattice. 

The enumeration of graphs contributing to Q(x) can be made easier by noticing that 
these graphs can be formed by combination of smaller graphs. We define an articulation 
point [11] as a point on y-axis such that the tree above is an allowed spiral tree in the 



quadrant above that part (note that these trees are defined in the upper quadrant and 
they never touch y = axis, except at (0,0)). For example, the solid squares represent the 
articulation points of the graph in Fig. 3, and Fig. 4 shows a spiral tree with no articulation 
point. Hence, these spiral trees can be seen as trees having y axis as a backbone on 
which spiral graphs are connected at different articulation points maintaining the spiral 
constraint. 

Let B(x) be the generating function of the quadrant spiral trees with no articulation 
points. Hence B(x) can be seen as sum of generating function of irreducible graphs with 
% sites along y-axis. We represent them by B^x) (see Fig. 5), then B(x) = J2i^i Bi(x). 
The full generating function in terms of B(x) would be 

Q(x) = x(l + B(x) + B\x) ) = l _ X B{x) (8) 

where B^x) are spiral graphs starting with i-sites along the y-axis. It is easy to see that 
Bi(x) = x, B 2 (x) = ^ and B 3 (x) = ^_ x _ x Ji^i_ x y One can write B^{x) with some effort 
but we do not have a general form for Bi(x) for any i. 

We restrict the graphs contributing to B^x) to be graphs such that they have % sites 
along the y axis and have at least one downward branch with % — 1 sites. This would not 
include structures like Fig. 4. We will represent the generating function of these graphs 
by Qi(x). Then we can represent B^x) in terms of two other generating functions, Vi(x) 
and Wi{x). We define Vi(x) as the generating function of spiral subgraphs starting with 
having % sites along y-axis. Wi(x) is the generating function of spiral subgraphs staring 
with i-sites along y-axis and ending with a downward branch with i — 1 sites (Fig. 6). 
Then, 

Bl(x) = Wi(x) + (9) 

x % 1 

Also, Vi(x) can be rewritten in terms of Wi(x) as 

V&) = xV^ix) + W^x) + ^MYizM (io) 

x l 1 

By expressing Qi(x) in terms of Bi(x) and Bi(x) in turn in terms of Wi(x), we can 
reduce the computational time. If W n is the number of graphs of size n contributing to 
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W(x) (W(x) = J2iliWi(x)), and Q n is the number of graphs of size n contributing to 
Qi(x), then W n grows more slowly than Q n . We enumerated W n and using them we could 
generate a 56 term series for Qi(x). The computation time for W n grows as (1.8) n , in 
contrast to (2.04) n for the Q n series. 

If we restrict the graphs contributing to Bi(x), Wi(x) and Vi(x) to the graphs having 
comb-like structure (by comb-like structure we mean graphs with one dimensional back- 
bone having vertical straight lines of arbitrary lengths), then it turns out that one can get 
the exact expression for these generating functions. We represent them by Vi(x),Wi(x) 
and Bi(x). It is easy to see that for comb like structures, 

W i (x)>W i (x) = - + - T ±1 + -, 7^ + 11 

1 — x 1 — xl — x 1 — x \1 — X J 

where K(x) = x 2 ^ l ~ 2 1 x j . Hence, 



Similarly, we get 



Wi(x) = x2%( \ X ] — (12) 

iV ; 1 - 2x + x 2 - x 3 + x l+1 V ; 



Vt{x) = l-2x + x 2 - X 3 + x^ (13) 



and hence 



x 2i (l-x) 2 

Bi(x) > B t (x) = ( 1 _ 2x + x2 _ x3 + xi ^ 1 _ 2x + x2 _ x3 + xi+1) (14) 

Substituting in Eq. 9 we get the generating function, Qi(x) for this subset of spiral 
trees in a quadrant. This generating function has a singularity at x c = 0.51662 which 
gives the growth constant A' of these trees to be 1.93565. Since this counts only a subset 
of all the spiral trees on a square lattice, this is a rigorous lower bound on \ sp i T ai for spiral 
trees on a square lattice. 

For the full Qi(x), we derived a 56 term series. If we assume, 

Q n ~ A>- ei (15) 
then we got estimates of Ai and 6\ to be 
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X 1 = 2.0449 ± 0.0001 



0i = 0.830 ±0.01 



(16) 
(17) 



1.2 Exact enumeration 

Since the number of configurations of a given cluster size is exponential in cluster size, 
the computational complexity of the algorithm for enumeration of all lattice animals or 
trees grows exponentially with the cluster size. For direct enumeration algorithms like 
Martin's algorithm [12], the time required to generate all the configurations of a given 
size grows as A n , where A is the growth constant and n is the cluster size and the memory 
requirement grows like a polynomial in cluster size. For lattice trees and animals, a finite 
lattice method [13] with an associated transfer matrix algorithm was used by Conway 
[14]. Conway generated a 25 term series for lattice animals using this algorithm. This 
series has recently been extended to 46 terms by Jensen [10] with a slight improvement in 
the algorithm. Both space and time requirements of this algorithm are found numerically 
to approximately grow as 1.4 n . The growth constant of lattice animals in contrast is 
approximately 4.06. Hence a considerable improvement in time is obtained by the transfer 
matrix algorithm at the cost of memory. 

The spiral constraint, cannot be easily implemented using the transfer matrix. Hence 
we have used Martin's algorithm for spiral trees, making use of the four-fold rotational 
symmetry of the lattice. Our series for number of trees and their average moment of 
inertia is given in Appendix. 1. 

Using this we generated a series of spiral trees on square lattice up to 37 terms (Ap- 
pendix 1). Earlier known series had only 25 terms. 

For analysing the series data we tried a four parameter sequential fit to the data of 
the form 



where 5 is an adjustable fixed parameter and B is a constant. We did a linear fit on the 
logarithm of Eq. 18 using A n , A n+ i, A n+2 and A n+3 to estimate values of B n , 5 n , \ n and 



A n = B\ n (n + 5) 



-9 



(18) 
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6 n . For spiral trees on square lattice we found a good convergence in successive values of 
A n and 6 n for 5 lying between 2.03 and 2.04. Fixing 5 = 2.0367 and B = 0.18124 we get a 
very good convergence of A n and 9 n for different values of n. These are given in Appendix 
2. From this we estimate 

A = 2.11433 ±0.00010 (19) 
9 = -1.3667 ±0.0010 (20) 

We have tried fits with non analytic corrections to scaling of the form , BX n (n+8) [1+ 
a/n A ], but we didn't get good convergence for A. Instead, B\ n {n + 5)~ e [\ — ae _/3n ] seems 
to fit much better with a m 0.32 and (3 0.35. 

For the radius of gyration data we used a sequential fit of the form 

log/ i>w = (2ui + l)ln(n + 5) + u + - V (21) 

(n + o) z 

where % stands for pi or ± as the case maybe and u and v are constants. 

For spiral trees in a plane Jj_,n would be zero and by symmetry the sum of square of 
x coordinate of all sites for all configurations of clusters of size n is symmetric with sum 
of squares of y-coordinate. Using Eq. 21 for sequential fit to our 35 term series we get 
a good convergence for 8 lying between —0.33 and —0.35. Fixing 5 = —0.338 we get the 
estimates of v vX to be 

2u lA = 1.3148 ± 0.0010 (22) 

These estimates are much more precise than the earlier estimates A = 2.1166 ± 0.001, 
= -1.307 ± 0.006 and 2v pl = 1.306 ± 0.010 using a 25 term series [5]. We can rule out 
the dimensional reduction conjecture with fair confidence. 

Above we presented our estimates using four parameter fits. Method of differential 
approximants has almost become a standard technique for such analysis [15]. In this case, 
the generating function has a divergent singularity at x c . We tried zeroth order differential 
approximants, they are listed in Table 1. We find a very poor convergence in values of 
x c and 9. Out of 70 approximants, 15 show spurious singularities (singularities with 
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\x c \ < 0.45). We have listed 20 values which showed best convergence. From these we get, 
A = 2.1142±0.002 and 6 = -1.39±0.02. Clearly the series is not very well behaved. This 
is reflected in the slow convergence of our series. Also Monte-Carlo generated random 
spiral trees of sizes 1000 (Fig. 2) suggest that the asymptotic behaviour of the series 
might set in rather late. Because of poor convergence of differential approximants, we 
have relied on parameter fits for series analysis in this paper. 

1.3 Monte- Carlo analysis 

With exact enumeration, we are restricted to clusters of size thirty seven in two dimen- 
sions. The main problem is with the extrapolation since the crossover sizes are likely to 
be large, as the total angle turned by the largest spiral arm about the origin for a spiral 
tree of size 1000 is about 2n only (Fig. 2). This indicates that the crossover value above 
which asymptotic behaviour sets in would be of order 10 3 . We tried to study larger spiral 
trees using MC methods. Monte-Carlo simulation of branched polymers is a challenging 
problem. Because of branching, most MC algorithms which are good for linear polymers 
show critical slowing down for branched polymers. For lattice trees there have been some 
studies using the cut and paste dynamic MC technique [9]. But with spiral constraint, 
algorithms involving large scale non local moves are not useful. We used an improved 
version of incomplete enumeration algorithm proposed recently by us [7]. Using it we 
could study spiral trees of sizes up to one thousand on a square lattice. 

Incomplete enumeration is a simple modification of exact enumeration algorithm and 
can be seen as a percolation process on the genealogical tree of the underlying enumeration 
problem. The optimal behaviour of the algorithm is achieved when we work around the 
percolation threshold of the genealogical tree. This algorithm falls in the class of stochastic 
growth algorithm like PERM [16]. We have shown in [7], that the asymptotic time to 
produce an independent sample of n sites for trees and animals grows as exp(an b ) with 
< b < 1 for this algorithm. Though the coefficient in front of stretched exponential can 
be made small by optimising the algorithm for small sizes. We will not give more details 
of the algorithm in this paper. These can be found in [7]. 

Fig. 2 shows pictures of some typical spiral trees of one thousand sites. Clearly, their 
structure is very different from lattice trees without the spiral constraint. Because of the 
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constraint they tend to branch much less. For spiral constraint, earlier numerical evidence 
suggest that unlike lattice trees and animals, spiral trees and animals do not lie in same 
universality class. The reason is that by allowing loops, the polymer can bend much more 
often and hence spiral animals would be more compact than the spiral trees. 

We studied spiral trees up to sizes 1000 using incomplete enumeration MC method. 
We made 10 7 MC runs. The moment of inertia tensor I p i yn as a function of n is plotted 
in Fig 7 and Fig 8. Assuming the asymptotic form to be such that 



In incomplete enumeration MC algorithm [7], each configuration of n sites is generated 
with equal probability P n which is just lT7=i where Pi is the probability with which 
an edge between level % and % + 1 on the genealogical tree of the problem is chosen. By 
keeping track of the average number of clusters of a given size generated in a given run, 
one can estimate the growth constant A and the critical exponent 9. But, the variance 
of the number of clusters increases as exp(n a ), < a < 1 for large n. Hence, instead we 
counted the number of descendants of each spiral tree generated. This approach has been 
used previously in [5, 17]. The mean number of descendants of a tree of size n gives a 
direct estimate of A n+1 /A n . We represent the mean number of descendants by M n . This 
is plotted in Fig. 9. A linear fit of the form A(l — 9/n) to this data gives A = 2.116 ± 0.01 
and 9 = —1.29 ± 0.02. For better estimates we assume 



With this we get the following estimates for n < 200 which are in agreement with the 
value obtained by extrapolating the exact series expansions. 



\og(I phn ) = logC + {2u p i + l)logn + — 
Using above written form, we get the estimate of v pi to be (Fig 7 and Fig 8) 



(23) 



2i/„, = 1.312 ±0.010 



(24) 




(25) 



A = 2.1145 ±0.0010 



9 = -1.364 ±0.010 



(26) 
(27) 



with 5 = 1.8. 
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2 Spiral trees on a cubic lattice 

In dimensions higher than two, the spiral constraint defined as the projection of path 
joining any site of the tree to the root in x — y plane contains no left turn can be employed 
in two ways. Bose et. al. [5] defined it such that for the projected path from origin to 
site on x — y plane only forward and right turns are allowed. But in dimensions higher 
than two, we can define another variation where trees as long as they do not violate the 
tree constraint and the projection on x — y plane is spiral, are allowed. We will call the 
spiral trees with only forward and right turns allowed as ST ± . 

If we allow for back-turns also, we would get different series because for example, 
Fig 10 shows one spiral tree of six sites which would not be a valid configuration if we 
consider only forward and right turns. We call the spiral trees with back-turns allowed as 
ST 2 . Naively, one would expect these two to belong to the same universality class. We 
generated the series till n — 17 on a cubic lattice using both definitions, however we find 
the two series behaving differently Series for both STx and ST 2 are given in Appendix 1. 

For STi, for A n the number of configurations, using Eq. 18 we find that the sequential 
fit shows a good convergence around 5 = 2.43. With 5 = 2.43 and B = 0.094, the values 
of A and 9 obtained are listed in Table 1. For u pi and u±, we used fitting form as given 
in Eq. 21, with 5 = —1.46 and 5 = —0.43 respectively. The sequential fits are given in 
Appendix 2 and estimates are listed in Table 2. 

Similarly, we obtained 17 term series for ST 2 . The sequential fits are given in Appendix 

2 and the values of A, 9, v v i and v± are listed in Table 2. 

The difference in value of A for STi and ST 2 is understandable as ST 2 has a greater 
number of configurations. More surprisingly, the critical exponents 9, v v i and uj_ within 
our error estimates are different in two models. In neither case, the conjectured dimen- 
sional reduction(Eq.l and 2) seems to be satisfied. 

3 Spiral trees in four dimensions 

On a hyper cubic lattice in four dimensions we generated a series till n — 13. We also 
correct mistakes in the earlier series reported for STi in [5]. The corrected series in given 
in the Appendix 1. We also obtained a 13 term series for ST 2 (see Appendix 1). The 
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estimates of A and critical exponents are listed in Table 2. 

We also performed Monte-Carlo simulations using incomplete enumeration algorithm 
for spiral trees up to size 50. Our estimates from MC simulations for STi are given in 
Table 3. 

Though we cannot rule out the possibility of 6 being zero in both series analysis and 
Monte-Carlo simulations, but it seems unlikely. 

4 Discussion 

Bose et. al. gave a plausible argument of curling up of the dimensions in the spiralling 
plane and had conjectured a dimensional reduction by four for spiral trees. Our nu- 
merical evidence as presented in this paper does not support the conjecture. The spiral 
constraint for trees seems to be very special. For example, the structure of spiral trees 
is very different from spiral animals with loops allowed [18]. Different implementation 
of the constraint in d > 2, seems to give different critical behaviour, suggesting different 
universality classes. A variety of self avoiding walks with different step restrictions rules 
on simple cubic lattice were studied in [19] using exact enumeration. Their analysis sug- 
gested same universality class for self avoiding walks with various restrictions (including 
the spiral constraint), as the unrestricted self avoiding walks. In contrast, our studies 
show different critical behaviour of spiral trees with different geometrical restrictions in 
three and four dimensions. 

We should note that for the large clusters of size 10 3 generated by Monte Carlo, the 
total angle turned by the largest spiral arm about the origin is about 2n. It is possible 
that the structure of spiral trees is such that this angle tends to infinity as n tends to 
infinity. In this case the crossover value above which asymptotic behaviour sets in would 
be expected to be of order 10 3 , and series analysis for smaller n may not give correct 
limiting behaviour. One indication that trees where spiral turns a lot are important is 
that the growth constant for spiral trees in a quadrant Q\{x) seems to be significantly 
smaller than for full spiral trees. 

For quadrant spiral trees on a square lattice, we obtained exact series up to sizes 56. 
There are very few such long series known for lattice models. The series gives a estimate 
of A = 2.044 for these quadrant spiral trees. This value is significantly smaller than for 
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the full spiral trees. 
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6 Appendix 1 

Exact series enumeration values in different dimensions 
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6.1 Two dimensions 



cluster size(n) 




(Ipl,n) 


1 


1 





2 


4 


1 


3 


14 


3.142857 


4 


40 


6.800000 


5 


105 


12.266667 


6 


268 


19.656716 


7 


674 


28.919881 


8 


1660 


40.159036 


9 


4021 


53.513056 


10 


9612 


69.074906 


11 


22734 


86.926014 


12 


53276 


107.140851 


13 


123916 


129.787372 


14 


286376 


154.926432 


15 


658100 


182.624835 


16 


1504900 


212.938547 


17 


3426464 


245.919131 


18 


7771444 


281.619675 


19 


17565064 


320.089299 


20 


39576360 


361.374917 


21 


88916877 


405.522760 


22 


199252252 


452.577078 


23 


445438310 


502.580546 


24 


993616344 


555.575100 


25 


2211923712 


611.601183 


26 


4914811468 


670.697934 


27 


10901498938 


732.903853 


28 


24141259980 


798.256392 


29 


53379537257 


866.791847 


30 


117861710196 


938.545859 


31 


259891311248 


1013.553288 


32 


572356464452 


1091.848086 


33 


1259008971656 


1173.463504 


34 


2766351037428 


1258.432171 


35 


6071954146120 


1346.786006 


36 


13314252070412 




37 


29167189621351 
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6.2 Three Dimensions ST X 



n 


A n 


(Ipl,n) 


(I±,n) 


1 


1 


0. 


0. 


2 


6 


0.66666 


0.333333 


3 


41 


1.85366 


1.07317 


4 


260 


3.63076 


2.27692 


5 


1568 


6.02296 


3.98214 


6 


9190 


9.06464 


6.19913 


7 


53090 


12.75954 


8.91987 


8 


303900 


17.09588 


12.1405 


9 


1727691 


22.0606 


15.8606 


10 


9767426 


27.6424 


20.0821 


11 


54966550 


33.8322 


24.8071 


12 


308138528 


40.6214 


30.0376 


13 


1721739000 


48.0022 


35.7754 


14 


9592901762 


55.9676 


42.0229 


15 


53314247488 


64.5112 


48.7822 


16 


295644339728 


73.6274 


56.0556 


17 


1636179620652 


83.3112 


63.8454 



6.3 Three Dimensions ST 2 



n 


A n 


(Ipl,n) 




1 


1 





0. 


2 


6 


0.666666 


0.333333 


3 


41 


1.85366 


1.07317 


4 


260 


3.63076 


2.27692 


5 


1576 


6.00762 


4.00761 


6 


9342 


9.00192 


6.30208 


7 


54890 


12.60084 


9.17041 


8 


320952 


16.7848 


12.6182 


9 


1869907 


21.5398 


16.651 


10 


10861750 


26.8572 


21.2772 


11 


62939998 


32.7312 


26.5047 


12 


364004296 


39.156 


32.3409 


13 


2101795408 


46.1276 


38.7927 


14 


12119643810 


53.6422 


45.8667 


15 


69805863940 


61.6968 


53.5693 


16 


401668665200 


70.2898 


61.9068 


17 


2309283532000 


79.4192 


70.8851 
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6.4 Four Dimensions ST\ 



n 


A n 


(Ipl,n) 


(I±,n) 


1 


1 


0. 





2 


8 


0.5 


0.5 


3 


80 


1.35 


1.5 


4 


800 


2.54 


3.030 


5 


7912 


4.05864 


5.10010 


6 


77656 


5.89816 


7.70862 


7 


759172 


8.04822 


10.84584 


8 


7403292 


10.49742 


14.50268 


9 


72073417 


13.23410 


18.67008 


10 


700774524 


16.24692 


23.34 


11 


6806914432 


19.52526 


28.5052 


12 


66064406668 


23.0592 


34.1596 


13 


640741734643 


26.8396 


40.2974 



6.5 Four dimensions ST2 



n 


A n 


(Ipl,n) 




1 


1 


0. 


0. 


2 


8 


0.5 


0.5 


3 


80 


1.35 


1.5 


4 


800 


2.54 


3.030 


5 


7960 


4.05226 


5.10754 


6 


79048 


5.87628 


7.74208 


7 


785748 


7.99822 


10.93174 


8 


7822676 


10.40506 


14.6724 


9 


78011513 


13.08484 


18.95778 


10 


779189988 


16.0274 


23.7816 


11 


7793589943 


19.22410 


29.1376 


12 


78049537766 


22.6676 


35.0206 


13 


782489000000 


26.3518 


41.4252 
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7 Appendix 2 



7.1 Two Dimensions Sequential Fit 



n 




n 




5 


2.078982187 


-1.4143616 


1.2918751 


6 


2.118727624 


-1.3598402 


1.3047319 


7 


2.117039314 


-1.3623751 


1.3092198 


8 


2.115352878 


-1.3651395 


1.3108492 


9 


2.114617151 


-1.3664433 


1.3117861 


10 


2.114771869 


-1.3661493 


1.312420 


11 


2.113813740 


-1.3680905 


1.3128895 


12 


2.113978775 


-1.3677359 


1.3132536 


13 


2.114183882 


-1.3672706 


1.3135423 


14 


2.114099443 


-1.3674721 


1.3137672 


15 


2.114103267 


-1.3674625 


1.3139586 


16 


2.114205656 


-1.3671946 


1.3141194 


17 


2.114223238 


-1.3671466 


1.3142505 


18 


2.114256310 


-1.3670527 


1.3143596 


19 


2.114279786 


-1.3669834 


1.3144497 


20 


2.114291286 


-1.3669483 


1.3145234 


21 


2.114301033 


-1.3669174 


1.3145839 


22 


2.114310834 


-1.3668854 


1.3146334 


23 


2.114311487 


-1.3668832 


1.3146734 


24 


2.114314464 


-1.3668728 


1.3147059 


25 


2.114318963 


-1.3668566 


1.3147321 


26 


2.114320428 


-1.3668513 


1.3147529 


27 


2.114321722 


-1.3668464 


1.3147694 


28 


2.114324605 


-1.3668351 


1.3147823 


29 


2.114326551 


-1.3668274 


1.3147921 


30 


2.114327932 


-1.3668217 


1.3147994 


31 


2.114329734 


-1.3668142 


1.3148047 


32 


2.114331349 


-1.3668072 


1.3148083 


33 


2.114332328 


-1.3668029 


1.3148104 


34 


2.114333055 


-1.3667997 


1.3148113 


35 


2.114333550 


-1.3667974 


1.3148113 


36 


2.114333553 


-1.3667974 




Es. Val. 


2.11433 ±0.0001 


-1.3667 ±0.001 


1.3148 ±0.001 
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7.2 Three Dimensions ST X 



n 




n 


2v p l,n 




5 


5.153269 


-1.019107 


0.847865 


1.128949 


6 


5.275382 


-0.810187 


0.865641 


1.098419 


7 


5.310873 


-0.743662 


0.871330 


1.083814 


8 


5.319667 


-0.725590 


0.874191 


1.077873 


9 


5.327658 


-0.707695 


0.875922 


1.073809 


10 


5.334141 


-0.691977 


0.876525 


1.070550 


11 


5.337903 


-0.682161 


0.876502 


1.068326 


12 


5.339533 


-0.677605 


0.876303 


1.067085 


13 


5.340111 


-0.675880 


0.876176 


1.066526 


14 


5.340282 


-0.675339 


0.876139 


1.066334 


15 


5.340255 


-0.675428 


0.876197 


1.066365 


Es. Val. 


5.340 ± 0.02 


-0.675 ±0.05 


0.876 ±0.05 


1.066 ±0.05 



7.3 Three Dimensions ST 2 



n 




On 




2^±,n 


5 


5.694072 


-0.143802 


1.054721 


1.200021 


6 


5.719085 


-0.123718 


1.013833 


1.171160 


7 


5.710159 


-0.132441 


0.989310 


1.153823 


8 


5.695471 


-0.149408 


0.977105 


1.147515 


9 


5.689143 


-0.157845 


0.969817 


1.144294 


10 


5.687350 


-0.160552 


0.963977 


1.141637 


11 


5.686110 


-0.162645 


0.959561 


1.139686 


12 


5.684763 


-0.165153 


0.956738 


1.138646 


13 


5.683809 


-0.167099 


0.955250 


1.138243 


14 


5.683473 


-0.167843 


0.954653 


1.138124 


15 


5.683632 


-0.167463 


0.954662 


1.138134 


Es. Val. 


5.683 ±0.02 


-0.167 ±0.05 


0.954 ±0.05 


1.138 ±0.05 
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Table Captions 



1. Table 1: Estimates of critical exponents and growth constant from differential ap- 
proximants. We looked at approximants for I > 9 and I — 3 < m < I + 3. We have 
tabulated here 20 values which showed best convergence. 

2. Table2:Estimates of critical exponents and growth constant from series analysis 
in three and four dimensions. Note that the value of 9 and v for rooted lattice 
animals/trees in 3d and Ad is known exactly(in 3d, 9 = v = 1/2 and in Ad, 9 = 5/6 
and v = 5/12 

3. Table3:Estimates of critical exponents and growth constants from Monte-Carlo sim- 
ulations in four dimensions. 
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Figure Captions 



1. Figurel:A rooted spiral tree of 15 sites on a square lattice. The root is the site 
enclosed in the square. At the root site the tree has freedom of choosing any of 
the four neighbouring sites. We count the spiral tree by number of sites and hence 
all bonds between two occupied sites is always assumed to be present. The site 
marked by X, if present will result in a loop for spiral site trees and hence will not 
be allowed. But it can be present in a spiral bond tree. 

2. Figure2: Randomly generated spiral trees of 1000 sites in 2-dimensions using incomplete- 
enumeration algorithm. 

3. Figure3:A simple counting problem of backbone with arbitrary long offshoots. Min- 
imum distance between two offshoots is 2 as else the tree constraint is violated. 
Solid squares represent the articulation points of the graph. 

4. Figure4:An example of an irreducible spiral graph with no articulation point. This 
is also an example of a graph not included in Qi(x) 

5. Figure5: Schematic figure of spiral trees contributing to Bi(x). Bi(x) is just a single 
vertex. 

6. Figure6: Example of graphs contributing to Vi(x) and Wi(x) respectively. 

7. Figure7:Plot of 4jrfe as a function of n for Monte-Carlo generated spiral trees on a 
square lattice. 

8. Figure8:Plot of I pljTl Vs n for Monte-Carlo generated spiral trees on a square lattice. 
The dotted line is a straight line with slope 2.312. 

9. Figure9:Monte-Carlo estimates of ratios of the number of configurations on a square 
lattice. The straight line gives a linear fit of the form A(l — 6/n) to the data. 

10. FigurelO:A spiral tree of six sites on a cubic lattice with a back-turn (drawn by a 
thicker line). This configuration will contribute to spiral trees ST 2 of six sites but 
not to STi. 
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[l,m] 


x c = 1/A 


n 
a 


[l,m] 


x c = 1/A 


n 
a 


[14, 13] 


n A 70000 V^d 

U.4 ( zooz5o 


— l.ooUoo 


[15, 18] 


n A 7QH71 /I /I 

OA tout 144 


— l.oyU to 


[14, 14J 


U.4 I ZyUo25 


— l.oooo4 


[14, 17] 


n /I 7on7ono 
U.4MU/oUo 


1 on 1 da 
— l.oylUo 


r T cr 1 Ql 

[15, 13J 


U.4 < zyU5io 


— i.oo4io 


[16, 17J 


U.4 < 60 1 006 


— l.oyzUy 


[16, 15] 


U.47294898 


-1.37U35 


[17,15] 


U.473U8675 


-1.39369 


[13, 15] 


U.47297513 


-1.37499 


[16,19] 


U.473U9U52 


-1.39421 


[16, 13] 


U.473U3UU7 


-1.384U9 


[18,15] 


U.4731U355 


-1.39686 


[13, 16] 


U.473U33U5 


-1.38436 


[17, 18] 


U.4731U9U6 


-1.39788 


[16,16] 


U.473U5593 


-1.388UU 


[15,16] 


U.47311UU1 


-1.39775 


[14, 15] 


U.473U5793 


-1.38863 


[18,18] 


U.47311U71 


-1.39822 


[15, 17] 


U.473U6712 


-1.39UU2 


[17, 19] 


U.47311U91 


-1.39826 



Table 1: Estimates of critical exponents and growth constant from differential approxi- 
mants. We looked at approximants for I > 9 and I — 3 < m < I + 3. We have tabulated 
here 2U values which showed best convergence. 





ST^d = 3) 


ST 2 (d = 3) 


STi(d = 4) 


ST 2 (d = A) 


A 


5.34U ± U.U2U 


5.683 ±U.U2U 


9.62±U.1U 


10.20 ±0.10 


e 


-U.675±U.U5U 


-U.167±U.U5U 


-0.11 ±0.10 


0.29 ±0.10 


Vpl 


U.44±U.U5 


U.477±U.U5 


0.34 ±0.05 


0.37 ±0.05 


v± 


U.54±U.U5 


U.69±U.U5 


0.44 ±0.05 


0.45 ±0.05 



Table 2: Estimates of critical exponents and growth constant from series analysis in three 
and four dimensions. Note that the value of 9 and v for rooted lattice animals/trees in 
3d and id is known exactly (in 3d, 9 — v — 1/2 and in Ad, 9 = 5/6 and v — 5/12 [6].) 





STi(d = 4) 


ST 2 (d = A) 


A 


9.60 ±0.1 


10.2 ±0.1 


9 


-0.13 ±0.1 


0.17 ±0.1 




0.33 ±0.02 


0.38 ±0.05 


v± 


0.451 ±0.020 


0.455 ±0.050 



Table 3: Estimates of critical exponents and growth constants from Monte-Carlo simula- 
tions in four dimensions. 
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Figure 1: A rooted spiral tree of 15 sites on a square lattice. The root is the site enclosed in 
the square. At the root site the tree has freedom of choosing any of the four neighbouring 
sites. We count the spiral tree by number of sites and hence all bonds between two 
occupied sites is always assumed to be present. The site marked by X, if present will 
result in a loop for spiral site trees and hence will not be allowed. But it can be present 
in a spiral bond tree. 
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Figure 2: Randomly generated spiral trees of 1000 sites in 2-dimensions using incomplete- 
enumeration algorithm 
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Figure 3: A simple counting problem of backbone with arbitrary long offshoots. Minimum 
distance between two offshoots is 2 as else the tree constraint is violated. Solid squares 
represent the articulation points of the graph. 




Figure 4: An example of an irreducible spiral graph with no articulation point. This is 
also an example of a graph not included in Q\{x) 
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Figure 5: Schematic figure of spiral trees contributing to Bi(x). Bi(x) is just a single 
vertex. 
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Figure 7: Plot of ^rm as a function of n for Monte-Carlo generated spiral trees on a 
square lattice. 
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Figure 8: Plot of I p i yn Vs n for Monte-Carlo generated spiral trees on a square lattice. 
The dotted line is a straight line with slope 2.312. 
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Figure 9: Monte-Carlo estimates of ratios of the number of configurations on a square 
lattice. The straight line gives a linear fit of the form A(l — 0/n) to the data. 




Figure 10: A spiral tree of six sites on a cubic lattice with a back-turn (drawn by a thicker 
line). This configuration will contribute to spiral trees ST 2 of six sites but not to STi. 
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